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ABSTRACT 

Engine-generator  units  for  use  in  standby  power  systems  are  reviewed. 
Reliability  models  including  state  transition  diagrams  are  developed 
for  sample  systems.  Methods  of  reliability  analysis  are  discussed  and 
ins  calculation  of  useful  reliability  measures  is  demonstrated.  A  computer 
program  for  the  solution  of  reliability  measures  is  included.  Example 
systems  studies  are  conducted  and  results  presented. 
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I.   INTRODUCTION 

The  inevitable  failure  of  electrical  power  networks  -  due  to  factors 
even  outside  the  system  manager's  control,  such  as  events  of  nature  and 
early  wear  out  of  components  -  dictates  the  necessity  for  stand-by  and 
emergency  power  generation  system  tc  support  critical  loads.  Many  factors 
can  influence  the  requirements  or  specifications  for  such  a  system,  but 
most  stem  from  an  evaluation  of  the  user's  needs.  Such  factors  include 
whether  or  not  life  threatening  or  hazard/safety  situations  are  involved, 
whether  plant  or  property  damage  could  occur,  whether  a  significant  loss 
in  revenue  from  operations  could  occur,  the  cost  of  back-up  equipment, 
and  the  reliability  of  the  normal  power  source.  This  paper  looks  at 
the  engine-driver  generator  system  as  an  emergency  or  standby  source 
and  investigates  a  method  for  modeling  the  integrated/standby  system 
for  the  purpose  of  obtaining  useful  reliability  information. 

Discussion  within  the  paper  begins  by  looking  at  typical  engine- 
generator  systems,  followed  by  a  presentation  of  techniques  for  evaluating 
reliability  measures  on  such  systems.  Reliability  models  and  state  transi- 
tion diagrams  are  then  created  for  example  engine-generator  configurations. 
Finally,  selected  parameters  within  the  models  are  varied  in  order  to 
observe  their  effect  upon  the  reliability  indices  and  the  results  are 
presented. 

II.  ENGINE-GENERATOR  CONFIGURATION 

Reference  [1]  discusses  many  arrangements  for  emergency  and  standby 
generation  systems.  Several  different  types  of  energy  conversion  hard- 
ware are  available  including  gas  and  diesel  engines,  gas  and  steam  turbines 
mechanical  stored  energy  systems,  and  battery  storage  and  uninterruptible 


power  supplies.  A  typical  engine-generator  configuration  is  shown  here 
in  Figure  1.  The  dotted  lines  show  how  the  system  can  be  expanded  to 
a  multi-engine  configuration.  There  are  many  advantages  to  installing 
several  units  of  lower  capacity  than  a  single  large  capacity  unit.  Be- 
sides the  increased  reliability,  which  shall  be  shown  later,  maintenance 
and  repairs  can  be  conducted  without  loss  of  adequate  backup,  additional 
less  critical  loads  can  be  supported  in  the  event  of  power  loss  while 
providing  the  necessary  degree  of  reliability  for  the  critical  loads, 
peak  load  showing  can  be  performed  with  extra  generating  capacity,  and 
future  expansion  is  facilitated  by  installation  of  similar  units.  Of 
course  the  advantages  must  be  balanced  against  the  costs  of  installing 
a  more  complex  system  as  well  as  the  length  of  time  the  generators  might 
be  exDected  to  operate  on  a  continuous  basis  -  larger  generators  be^g 
more  suited  to  longer  operating  periods. 

Operation  of  the  system  in  Figure  1  is  assumed  to  be  fully  auto- 
matic. When  the  normal  source  fails,  all  engines  in  the  system  are  automatically 
started.  The  first  engine  to  reach  operating  voltage  and  frequency 
causes  the  emergency  load  to  transfer  to  the  generator.  Other  engines 
are  paralled  and  switched  on-line  to  share  the  load.   If  sufficient  capa- 
city is  available,  other  less  critical  loads  can  then  also  be  transferred 
to  the  generators.  In  the  event  a  generator  should  fail,  appropriate 
load  shedding  automatically  occurs.  Regularaly  scheduled  inspections 
and  preventive  maintenance  are  assumed  to  take  place. 

An  important  matter  to  bear  in  mind  when  considering  the  application 
of  a  standby/emergency  system  is  to  ensure  loads  are  compatible  with 
the  monetary  interruptions  and  transients  associated  with  the  transfer 
between  power  sources.  Lighting  and  other  resistive  loads  are  highly 


tolerant  of  breaks  and  large  transients,  whereas  transients  lasting  on 
the  order  of  micro-seconds  could  be  disastrous  to  data  processing  and 
comunications  hardware -important  information  or  critical  program  execu- 
tion and  memory  could  be  lost.  Special  consideration  must  also  be  given 
to  highly  inductive  loads  and  associated  motor  control  circuitry. 

COMPONENT  FAILURE  MODES  AND  DATA  COLLECTION 

Before  delving  into  reliability  models  and  measures,  a  brief  dis- 
cussion on  statistical  data  collection  is  in  order.  The  collection  of 
information  concerning  equipment  failures  -  the  time  of  failure,  time 
to  repair,  severity,  and  cause  of  failure,  etc.  -  is  the  foundation  upon 
which  all  quantitative  reliability  analysis  rests.  The  statistics  generated 
from  this  data  in  the  form  of  mean  durations,  failure  and  repair  rates, 
probabilities,  etc.  are  the  basic  quantities  used  in  the  modeling  process. 
Reference  [2]  presents  several  discussions  which  highlight  the  importance 

and  fledgling  data  reports  and  analysis  systems.  The  engineer  in  the 
field  plays  a  vital  role  in  supporting  these  systems  by  ensuring  accurate 
and  timely  maintenance  and  repair  data  are  reported. 

When  adequate  data  bases  are  not  available,  more  subjective  techniques 
such  as  failure  modes,  effects,  and  criticality  analysis  must  be  employed. 
Significant  system  components  are  enumerated  and  the  modes  or  manner 
in  which  the  component's  failure  affects  the  system  are  listed.  A  sub- 
jective ordering  of  importance  of  this  list  in  accordance  with  some  appli- 
cable criteria  (life-safety,  endurance,  cost,  etc.)  can  be  followed  by 
non-parametric  statistical  techniques  to  provide  rudimentary  quantitative 
data.  A  listing  of  component  failure  modes  can  also  be  helpful  in  creating 


the  framework  for  a  data  collection  and  analysis  system  capable  of 
providing  accurate  statistics  for  more  complex  quantitative  reliability 
models;  an  example  of  this  could  be  a  model  which  distinguishes  the  dif- 
ference between  long  and  short  term  failures,  with  respective  frequencies, 
yielding  reliability  indices  more  accurate  than  results  from  a  lumped 
average  model . 

IV.  RELIABILITY  EVALUATION 

A.  BASIC  TERMS  AND  INTRODUCTION 

To  effectively  perform  reliability  studies,  it  is  helpful  for  the 
analyst  to  possess  a  working  knowledge  of  probability  theory  and  portions 
of  statistical  methodology.  Reference  [3]  provides  a  useful  introduction 
to  these  topics  and  employs  many  examples  directly  applicable  to  power 
system  studies.  Given  this  background,  a  few  additional  terms  and  defi- 
nitions are  required  before  an  analysis  can  be  undertaken.  Reference 
|_tj  nas  ccmpi  led  many  Oi  tne  commonly  usee  re<iau"i  1 1  lv  cerms;  pertinent 
trms  and  indices  for  this  paper  are  presented  here. 

Failure  Rate  -  mean  number  of  failures  per  unit  exposure  time  of 
a  component. 

Unavailabil ity  -  the  steady  state  probability  that  a  system  will 
be  in  a  failed  state  or  out  of  service  for  scheduled  maintenance,  etc. 
This  is  the  long  run  fraction  of  time  spent  in  the  failed  state. 

Frequency  of  System  Failure  -  mean  number  of  system  failures  per 
unit  time. 

Mean  Down  Time  -  expected  or  long-term  average  duration  of  a  single 
failure  unit. 


Mean  Cycle  Time  -  expected  time  between  successive  failures. 

Mean  Up  Time  -  expected  time  of  stay  in  the  up  state  in  one  cycle. 

B.  METHODS  OF  SOLUTION 

There  are  several  well-established  techniques  for  the  solution 
of  system  reliability  indices.  Chief  amongst  these  are:  network  tech- 
niques including  network  reduction  and  the  minimal  cut  set  approach, 
useful  when  a  system  can  be  represented  as  a  collection  of  series  and/or 
parallel  components;  decomposition  by  conditional  probability,  a  manner 
of  reducing  a  complicated  system  into  simple  subsystems  whose  reliability 
indices  can  be  figured  separately  and  combined  with  the  remaining  sub- 
systems; and  the  state  space  approach,  in  which  a  system  can  be  represented 
by  a  collection  of  states,  each  of  which  completely  describes  the  conditions 
under  which  the  system  is  operating,  e.g.  components  up,  components  down, 
the  given  environment,  etc.  The  state  space  approach  is  considered  ideal 
for  the  solution  of  the  standby/emergency  generator  problem  since  the 
system  states  can  be  quickly  enumerated  and  interdependences  associated 
with  component  failures  or  environmental  transitions  can  be  readily  accomo- 
dated. 

C.  STATE  SPACE  APPROACH 

Use  of  the  state  space  approach  for  the  solution  of  system  relia- 
bility indices  requires  four  primary  steps  as  developed  in  Reference  [5] 
and  shown  as  follows: 

1.  Identify  all  system  states  -  typically,  components  can  be 
operating  or  failed,  or  a  particular  environmental  condition 
can  be  in  effect. 


2.  Determine  interestate  transition  rates  -  the  mean  number  of 
occurrences,  per  unit  time,  in  which  the  system  undergoes  the 
change  or  transition  represented  between  state  i  and  state  j . 
This  figure  is  a  key  statistic  obtained  from  the  field  data 
collection  process.  An  important  assumption  made  here,  which 
is  r.ecassary   for  the  subsequent  analysis,  is  that  the  length 
of  stay  in  a  given  state,  considered  to  be  a  random  variable, 
possesses  an  exponential  probability  distribution.  The  prime 
consequence  of  which  is  that  the  interstate  transition  rates 
are  constant. 

3.  Computer  state  probabilities  -  the  probability  of  the  system 
being  in  state  i.  Here,  the  concept  of  statistical  independence 
comes  into  play.   If  a  component  is  allowed  to  undergo  any  one 
of  its  possible  transitions  out  of  a  given  state  without  being 
influenced  by  the  condition  of  the  remaining  components  within 
the  system,  or  the  prevailing  environment,  then  it  is  considered 
independent  of  the  other  components  and  the  environment.  When 
all  components  are  independent,  then  the  determination  of  the 
state  probability  becomes  a  simple  matter  of  computing  the  product 
of  the  component  probabilities  for  their  respective  conditions. 
If,  however,  there  is  statistical  dependence  amongst  the  compo- 
nents, e.g.  no  other  components  may  fail  after  one  has  failed, 

or  upon  the  environment,  then  a  set  of  simultaneous  linear  equations 
must  be  solved.  This  set  of  equations  stems  from  the  theory 
of  continuous  Markov  processes  (see  Reference  [3]);  this  is 
where  the  constant  transition  rate  assumption  is  necessary. 


Reference  [6]  presents  an  efficient  matrix  algorithm  for  the 

computation  of  state  probabilities  for  this  situation: 

BP  =  C  (1) 

B  =  matrix  obtained  from  A  by  replacing  all   elements 
of  an  arbitrarily  selected  row  !<  by  1; 

A  =  matrix  of  transition  rates   such  that  element 

aij  =  Xji  and  aii  =  -Zaji     . 

j,jfi    ' 

Aij  =  constant  transition  rate  from  state  i  to  j; 

P  =  column  vector  of  steady  state  system  state  probabilities; 

C  =  column  vector  with  Kth  element  equal  to  1  and  other 
elements  set  to  zero. 

Calculate  reliability  indices  -  Reference  [6]  provides  convenient 

formulas  for  calculation  of  the  most  commonly  used  reliability 

measures : 

(a)  system  unavailability  -  the  summation  of  probabilities 

of  the  state  representing  system  failure, 

Pf  =  E  pi  (2) 

T   ieF 

where  Px  =  probability  of  system  failure, 
F  =  subset  of  failed  states, 
pi  =  probability  for  state  i; 

(b)  frequency  of  system  failure  - 

ff     £    pi   £   Aij  ,  (3) 

ie(S-F)    jeF 

S   pi    Z    Aij  ,  (4) 

ieF     je(S-F) 


where  ff  =  frequency  of  system  failure, 
S  =  system  state  space; 

(c)  mean  cycle  time  - 

Tf  =  l/ff     ;     ff  -    (fi)  (5) 

(d)  mean  down  time  -  the  expected  time  of  stay  in  F  in  one-  cycle 
is 

i  [7  t  ;  (6; 

(e)  mean  up  time  -  the  expected  time  of  stay  in  up  states  (S-F)  is 

du  =  Tf  -  df  .  (7) 

It  should  be  noted  that  the  above  procedure  for  calculating  the 
probabilities  and  indices  was  based  on  steady-state  analysis.  Time 
dependent  solution  in  the  form  of  coupled  differential  equations  is 
also  possible,  albeit  more  complex,  should  transient  analysis  be  required. 
As  may  have  been  observed,  and  experience  bears  out,  the  state  space 
appraoch  can  become  an  unwieldly  process  for  large  interdependent  systems. 
It  proves,  though,  .to  be  an  ideal  technique  for  the  standby/emergency 
generator  problem. 

V.   RELIABILITY  MODELS  FOR  THE  EMERGENCY/STANDBY  GENERATOR  SYSTEM 

The  procedure  for  evaluating  reliability  measures  will  be  demon- 
strated through  the  use  of  two  separate  models.  Reference  [6]  developed 
a  simplified  reliability  model  for  a  combined  normal  and  standby/emergency 
power  source.  The  first  model  presented  here  is  an  expansion  on  that 
system.  The  second  model  used  in  this  paper  is  developed  for  a  dual 
generator  standby  configuration.   Both  of  these  models  are  based  on 


the  engine  generator  configuration  shown  in  Figure  1.  In  addition, 
these  models  shall  also  account  for  a  two-state  weather  environment. 
Figures  2  and  3  are  the  state  transition  diagrams  for  these  systems. 

In  the  diagrams,  all  the  possible  system  states  have  been  enumerated, 
E=:h  state  represents  the  following  conditions: 

a)  Normal  source  - 

NG  =  Normal  Good 
NF  =  Normal  Failed 

b)  Standby  generator(s)  - 

1.  single  generator  case 

SG  =  Standby  Good 
SF  =  Standby  Failed 

Distinction  is  also  made  as  to  whether  the  standby  failed 
on  start-up  or  subsequently  while  running  under  load. 

2.  dual  generator  case 

G  =  generaxor  Is)   Guuu  ('either  one'or  bothj 
FR  =  generator(s)  Failed  while  Running 
FS  =  generator(s)  Failed  on  Start-up 

c)  Environment  or  weather  state  - 

Environment  #1  =  Normal  Weather 

Environment  #2  =  Inclement  (or  Stormy)  Weather 

The  interstate  rates  represented  on  the  diagrams  are  as  follows: 

Ap,  yp  =  Failure  and  repair  rates  for  normal  source 

As    =  Failure  of  standby  generator 

yss   =  repair  rate  of  standby  generator  for  start-up  failure 
ysr   =  repair  rate  of  standby  generator  for  running  failure 
co-, 2  w21  =  transition  rates  between  environments 

Unprimed  rates  denote  values  for  normal  weather  state,  primed  rates 
denote  values  for  inclement  weather  state.  With  the  distinction  of 
start  and  running  failures,  the  probability  of  starting  the  standby 
source  when  needed  must  be  accounted  for  by  the  coefficient  a.  The 
probability  of  a  failure  on  start-up  is  the  complementary  coefficient: 
a  =  1  -  a  . 
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Some  important  assumptions  must  be  mentioned  regarding  the  state 
transition  diagrams: 

1)  The  probability  of  two  failures  and/or  repairs  occuring  simul- 
taneously is  considered  negligible.  Hence,  each  interstate 
transition  represents  only  a  single  event. 

2)  All  repairs  are  considered  to  be  independent,  i.e.  if  the  normal 
and/or  one  or  more  generators  have  failed,  repairs  progress 
simultaneously  and  respond  according  to  their  respective 
transition  rates. 

3)  Generator  fuel  supplies  are  considered  inexhaustible,  thus 
lack  of  fuel  is  not  counted  as  failure. 

Treatment  of  Two-State  Environments 

References  [5]  and  [7]  discuss  a  method  for  evaluating  reliability 

measures  under  fluctuating  environments.  The  key  to  understanding  this 

technique  is  in  recognizing  the  interdepend es  created  between  the  duration 

of  stay  in  an  environment  an~  the  component  transition  rates  within 

that  environment.  The  following  relations  can  be  used  to  obtain  values 

for  transition  rates  within  a  given  environment. 

AAV  =  A[N/(N+S)]  +  A'  [S/(N+S)]  (8) 

A'S/(AN  +  A'S)  =  x  (90 

where 

AAV  =  average  failure  rate 

A  =  failure  rate  within  environment  #1 
A'  =  failure  rate  within  environment  #2 
N  =  length  of  stay  in  environment  #1 
S  =  length  of  stay  in  environment  #2 
X  =  fraction  of  total  number  of  failures  occurring  in 
environment  #2 
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Equation  (3)  relates  the  average  or  equivalent  failure  rate  (typically 
available  from  data  collection  activities)  to  the  respective  failure 
rates  within  each  environment,  each  being  multiplied  by  the  fraction 
of  time  the  system  stays  within  that  environment.  Equation  (9)  is  just 
a  manner  of  representing  the  fraction  of  total  failures  which  occur 
in  environment  #2.  The  transition  rates  between  environments  are  found 

by  inverting  the  lengths  os  stay  within  the  environments: 

i 
u1?  ~tt=   Transition  rate  trom  environment  one  to  two; 

oop,  =  -f  =  Transition  rate  from  environment  two  to  one. 

This  assumes  that  the  durations  of  the  environments  are  exponentially 
distributed,  thus  allowing  constant  transition  rates.  For  the  standby 
Generation  system  models  develcpeH  in  this  .sac-tion^it   has  also  been 
assumed  that  the  repair  rates  for  all  components  are  unaffected  by  the 
fluctuating  environment. 

It  should  be  mentioned  that  the  above  treatment  of  a  two-state 
environment  can  be  used  for  many  different  types  of  environments.  Besides 
weather  conditions,  other  environments  such  as  variable  corrosive  conditions, 
dust  conditions,  or  electromagnetic  field  conditions  could  also  be  con- 
sidered. 

VI.  EXAMPLE  SYSTEM  STUDIES 

For  the  transition  diagrams  of  Figures  (2)  and  (3),  reliability 
indices  will  be  obtained  by  computer  solution  in  a  manner  which  allows 
variation  of  selected  parameters  (rates  or  probabilities)  in  order  to 
observe  their  effects  upon  the  system. 
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A.  SAMPLE  DATA 

The  sample  data  used  in  Reference  [6]  for  the  engine-generator 
system  will  be  used  in  these  studies.  For  the  normal  source,  the  failure 
rate  and  mean  down  times  are 

AP(AV)  =  0.537 /year 
MDT^  =  mean  down  time 

u 

=5.66  hours 

i 
Note;     The  reciprocal  of  MDTn  is  unJ  thus  un  =  u^      =  1548.76325/year. 

For  each  engine  generator 
K(j\\j\   =  0.00536/year 
MDTS  =  478  hours 

ySR  =  18.33891/year 
From  reference  [5],  typical   values   for  the  mean  durations  of  each  environ- 
ment are 

r-«QPri       niivafinn      q  »-» \  /  i  v*r\  v-\rrtr\  v*>  ■»-       ■#  1  M      —      0("^(^       i^-  — .  .  .  _-  c- 

IIICUI I        U  lil     C*  O    i  W  i  i        Cil  *     i   i    vilmCiii.        i,    J-    5        M  L.KJ  \J        llUUi    J 

mean  duration  environment  #2,  S  =  1.5  hours 
Thus, 

uj12  =  43.830/year 

co21  =  5844.0/year 

B.  PARAMETER  STUDIES  AND  PROGRAMMED  SOLUTION 

The  program  used  for  solution  of  the  sample  systems  is  provided 
in  the  appendix.  A  series  of  four  studies  were  conducted  on  both  the 
single  engine  and  dual  engine  . 
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The  series  of  tests  were  repeated  for  the  generator  systems  in 
three  separate  environmental  configurations.  The  three  environmental 
configurations  were:  a)  one  state  environment  -  no  fluctuating  weather 
conditions  considered;  b)  two-state  weather  environment  with  the  fraction 
of  failures  in  the  stcrmy  state,  x  =  0.20;  and  c)  two-state  weather 
environment,  ;<  =  0.80.  The  series  of  four  studies  included  the  following: 

1)  Fix  u'  equal  to  ;.-R  and  vary  a  from  0.02  to  0.10; 

2)  Fix  a  equal  to  0.01  and  vary  y^c  from  yCD  to  2y<;nJ 

3)  Repeat  (2)  with  a  =  0.05; 

4)  Repeat  (2)  with  a  =  0.10. 

For  the  two-state  environment  configurations,  the  average  failure  rates 
were  broken  down  as  follows: 

=  0.20  =  0.80 


Ap  =  0.43282/yr 

XPI  =  14.4274/vr 
P 

Xs  =  0.00432/yr 

X$l  =  0.14401/yr 


Ap  =  0.10821/yr 
AD!    =  57. 70960 /yr 

A$  =  0.00108/yr 
Xs,  =  0.57602/yr 


C.   RESULTS 

Results  of  the  system  studies  are  presented  in  Tables  1-8.      Each 
table  contains  the  results  of  one  of  the  generator  systems   in  all 
three  environmental   situations.     The  following  observations  and  conclu- 
sions can  be  drawn  from  these  results. 

1.  In  the  case  of  the  single  generator/single  state  system,  with 
^SS  =  ^SR  ^e  resu^s  are  identical  to  those  obtained  for  the 
simplified  engine-generator  model   of  Reference  [6]. 
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2.  The  most  significant  improvements  seen  in  the  reliability 
indices  come  as  a  result  of  introducing  the  redundancy  of  the 
dual  generator  system.  Improvement  over  the  single  generator 
system  appears,  though,  to  be  dependent  upon  the  value  of 

a.  As  a  increases,  the  degree  of  improvement  becomes  less. 
At  a   =  0.01  there  is  a  factor  of  100  improvement  over  the  single 
engine  system,  while  at  a  =0.05  a  factor  of  only  20  exists 
and  at  a  =  0.10  a  factor  of  only  10  exists. 

3.  In  all  situtations,  the  indices  are  found  to  suffer  as  &   is 
increased.  For  the  single  generator  system,  the  increase  in 
expected  down  time  per  year  and  failure  frequency  over  the 
range  of  a  is  approximately  an  order  of  magnitude,  while  the-' 
indices  vary  two  order  ot  magnitude  for  the  dual  generator 
system.  The  variation  for  the  dual  generator  system  appears 
to  be  the  square  of  the  variation  for  the  single  generator 
system. 

4.  The  variation  of  y<.<-  over  its  range  resulted  in  only  small 
effect  on  the  reliability  indices  for  all  configurations. 
Although  slight,  the  change  was  noted  always  as  an  improvement 
in  the  indices  as  the  value  of  y~~  increased. 

5.  Introduction  of  the  weather  states  also  has  only  a  limited 
effect  on  the  indices.  There  is  mixed  results  as  to  whether 
improvement  on  degradation  of  the  indices  occurs.  It  appears 
that  a   plays  a  role  in  defining  a  threshhold  at  which  the 
two-state  environments  show  improvement  of  the  indices  over 
the  single  state  environment.  For  the  single  engine  system, 
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improvement  is  noted  for  a   >  0.01.  For  the  dual  generator 
system,  improvement  occurs  for  a  >  0.05.  The  change  in  the 
value  of  x  does  not  alter  these  trends,  but  does  increase  the 
magnitude  of  the  difference. 

CONCLUSION 

This  paper  has  discussed  the  engine-generator  in  the  role  of  a 
standby/emergency  power  system.  The  state  space  approach  for  evaluating 
important  reliability  measures  can  be  quickly  and  efficiently  applied 
to  varied  standby/emergency  configurations.  The  results  of  such  a  quan- 
titative evaluation  can  assist  in  making  cost  vs.  benefit  decisions 
for  an  installation  as  well  as  making  comparisons  of  generator  systems. 
The  example  studies  conducted  in  the  paper  demonstrated  the  improved 
reliability  measures  obtained  by  introducing  redundancy  into  a  standby 
configuration,  but  another  important  aspect  to  consider  is  that  of  the 
start-up  failure  probability.  Results  showed  that  a  single  engine  system 
with  high  start-up  probability  stands  to  be  as  equally  reliable  as  a 
dual  generator  configuration  with  lower  start-up  probability.  Consideration 
of  environmental  conditions  can  also  have  an  effect  on  reliability  measures 
but  on  a  far  more  restricted  scale. 
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TABLE 


I.  SOLUTION  OF  RELIABILITY   INDICES   FOF    SINGLE   GENERATOR 

mu(ss)   FIXED   8    18.33891,   ALPHA  VARIED 


CNE  STATE  ENVIRONMENT. 

TEST 

EXP  DOWN 

FAILURE 

MEAN 

PARAMETER 

TIME,  YEAR 

"E2LE-;:- 

CSXLE  TIME 

ALPHA 

HOURS  ) 

EAR 

EARS  )    i 

;■  j  '  0000 

0.030893 

0.065524 

181-03709 

8.820068 

0.061  "c3 

8.0  i  '  042 

90.5663? 

0.836068 

0.692520 

0.016558 

60 

39i53 

0.04-0000 

0.123455 

0.e22070 

45 

31030 

0.056000 

0.1542"-. 

8  02_573 

36 

25334 

0.066906 

0.185070 

0  033085 

30 

22519 

8.078898 

0.21  5S50 

0.03S5S8 

25 

91500 

0.086060 

0.246613 

6.044087 

22 

68232 

0.090600 

0.277359 

0.049584 

20 

16797 

0.106000 

0.308086 

6.055677 

18 

15648 

MEAN  MEAN 

DOWN  TIME  UP   TIME 

HOURS   )  (   YEARS    ) 


5.594 

181.63645 

5.59A 

96.55375 

5.53a 

60.39339 

5.594 

i5. 30966 

5.594 

36.25870 

5.594 

30.22456 

5.594 

25.91437 

5.594 

22.68168 

5.594 

20.16734 

5.594 

18.15585 

TWO-STATE   ENVIRONMENT: 


TEST 

EXP.  DOWN 

FAILURE 

MEAN 

PARAMETER 

TIME/YEAR 

FREOUENCY 

CYCLE  TIME 

ALPHA 

(  HOURS  ) 

(1/YEAR) 

YEARS  )    ( 

0  010000 

0.030901 

0.005524 

181.02392 

0 

026000 

0.061762 

0 

01 1041 

96.57022 

0 

030000 

0.092605 

0 

016555 

60.40478 

0 

040000 

0.123430 

0 

022066 

45.31938 

0 

0500O0 

0.154237 

0 

027573 

35.26730 

0 

060000 

0.185027 

0 

033077 

30.23222 

0 

070600 

0.215793 

0 

0365-3 

25.92127 

0 

080000 

0.246552 

0 

044076 

22.68796 

0 

090000 

0.277233 

0 

049571 

20.17312 

c 

1CC0C0 

o 

055062 

15.16120 

MEAN  MEAN 

DOWN   TIME  UP   TIME 

(   HOURS   )  (  YEARS   ) 


5.594 

181.02328 

5.594 

90.56958 

5.594 

60.40414 

5.594 

45.31875 

5.594 

36.26666 

5.594 

30.23158 

5.594 

25.92063 

5.594 

22.68733 

5.594 

20.17248 

Cj.oyt 

16.16656 

TWO-STATE   ENVIRONMENT:   x 


TEST 

EXP.  DOWN 

FAILURE 

MEAN 

MEAN 

MEAN 

PARAMETER 

TIME/YEAR 

FREOUENCY 

CYCLE  TIME 

DOWN  TIME 

UP  TIME 

ALPHA 

(  HOURS  ) 

(1/YEAR) 

YEARS  )    ( 

HOURS  )    ( 

YEARS  ) 

0.018000 

0.030933 

0.005531 

180.80253 

5.594 

180.88198 

0.020000 

0.061660 

8.81  1023 

98.7281  1 

5.594 

90.71947 

0.030000 

0.092362 

0.016512 

60.56347 

5.594 

68.56283 

0.040000 

8.123846 

0.021997 

45.46082 

5.594 

45.46818 

0.050800 

0.15371  1 

0.027479 

36.39142 

5.594 

36.39878 

0.060000 

0.184358 

8.032958 

30.34191 

5.594 

30.34127 

0  670000 

0.214936 

6  038433 

26.01923 

5.594 

26.01859 

0  088006 

6.245595 

0.043905 

22.77634 

5.594 

22.77571 

0.096O00 

0.276186 

0.849374 

20.25353 

5.594 

20.25295 

0.100800 

0.306759 

8.054839 

18.23505 

5.594 

13.23441 

TABLE   2. 


SOLUTION   OF  RELIABILITY   INDICES   FOR   DUAL  GENERATOR 
mu(ss)   =    18.32391.   aipna   VARIED 


ONE   STATE   EN  /IRONMENT 


"5- 

:  dan 

FAILURE 

MEAN 

MEAN 

MEAN 

PARAMETER 

T  '-r   EAR 

-redle:::- 

CYCLE  TIME 

DOWf' 

1  TIME 

UP  TIME 

alpha 

HOURS  ) 

■     rEAR 

f£  AR  S 

HOURS 

)    ( 

YEARS  ) 

0.010008 

0.00031  S 

0.000053 

17364.452  '  3 

C  =-c 

17364.45150 

e  020090 

e.0012*: 

0.00023  1 

4346.  :4561 

5.5:3 

4346.24498 

0.020090 

0.002860 

0.0005  '  "" 

'  933. 50802 

5.529 

1953.50739 

8  040000 

0.005079 

0.060?  '  9 

1088.57170 

5.529 

1088.57107 

0.050000 

0.00-;:- 

0.001  434 

697.291  71 

5.529 

697.29108 

0.060000 

0.01  14.08 

0.002063 

484. 64456 

5.529 

434  64393 

0.070000 

0.015515 

0.002806 

356.36663 

5.529 

356.36600 

0.030000 

0.020248 

0.803662 

273.07213 

5.529 

273.07150 

0  O900O0 

0.025605 

0.004631 

215.94052 

5.529 

215.93989 

0.100000 

0.031584 

0.005712 

175.05676 

5.529 

175.05613 

TWO-STATE   ENVIRONMENT: 


0.20 


TEST 

EXP.  DOWN 

FAILURE 

MEAN 

PARAMETER 

TIME/YEAR 

FREQUENCY 

CYCLE  TIME 

alpha 

[  HOURS  ) 

(1/YEAR)    (  YEARS  )    ( 

0.010000 

0.000319 

0.000058 

17354.35471 

0.9200O0 

8.001272 

0.800230 

4345.27835 

0.03U000 

0.002860 

0.000517 

1933.31334 

0.940000 

8.095079 

0.000919 

1088.53021 

0.050000 

8.097929 

8.0O1434 

6S7. 29223 

0.060000 

3.91  1408 

0.002063 

434.65798 

0.070000 

0.015514 

9.002306 

356.28269 

0.038000 

6.9202A6 

0.093562 

273.09956 

0.390000 

0.925603 

9.904631 

215.95713 

P  'PC-PC* 

C3 ~ 1 5S2 

CCC5-12 

175.07216 

TWO-STATE  ENVIRONMENT: 

x  =  0.80 

MEAN  MEAN 

DOWN   TIME         UP  TIME 
HOURS   )         (  YEiRS   ) 


5.529 
5.529 
5  529 
5.529 
5.52S 
5.529 
5.529 
5.529 
5.529 


17354.35408 
4345.27772 
1933.31271 
1088.52958 
697.29160 
434.65735 
356.33306 
273.08893 
215.956*9 
175.37153 


TEST  EXP.   DOWN  FAILURE  WEAN  MEAN  MEAN 

PARAMETER         TIME/YEAR  FREOUENCY       CYCLE  TIME  DOWN   TIME  UP  TIME 

olpho  (   HOURS   )  (1/YEAR)  (   YEARS   )         (   HOURS   )  (   YEARS   ) 


0.810000 

0.920000 
0.030000 
0.048000 
0.050080 
0.0600OO 
0.878000 

e  080000 

0.090000 
0.100000 


0.00 
0.00 
0.00 
0.00 
0.01 
0.01 

0.02 

0.02 
0.03 


1277 
2355 
5033 
7939 

1404 

5504 

0227 
5573 
1539 


0.000058   17193.77553 
0.000231    4329.65876 


0.000518 
0.000919 
0.001434 
0.002063 
0.002804 
0.003658 
0.004625 
0.085704 


1930.09830 
1087.80185 
697.25629 
484.84096 
356.53270 
273.34959 
216.20773 
175.30623 


5.529 

17193.77500 

5.529 

4329.65813 

5.529 

1930.09767 

5.529 

1087.80122 

5.529 

697  25566 

5.529 

484.84033 

5.529 

356.63207 

5.529 

275.34396 

5.529 

216.20710 

5.529 

175.30560 

TABLE 


SOLUTION  CF  RELIABILITY  INDICES   rOR   SINGLE 
clcho   =   0.01  mu(ss)  voried 


JENERATOR 


DNE   ST-"E   ENVIRONMEh 


TEST 

?aRAME~E 
muiss'i 


'  ■-     DC 

EAR 


WLURE 
FREOUENCY 
EAR] 


MEAN 

CYCLE  TIME 
fEARS    ■ 


MEAN  MEAN 

DOWN  TIME         UP  TIME 
HOURS   )  (  YEARS   ) 


■5.333?'  a 

7  280 1 
Z2.066692 
23.840533 
25.67447* 
27.50836£ 
29.342251 
31. 1 76- 4~ 
33.010038 
34.843929 
36.677820 


"  ?39B 

330783 

3366S2 

9  330596 

-  30507 

30431 

o. 330:3^ 

30292 
0.030229 
0.030169 

0.0301  1  1 


S5524 

15510 

•  35498 

0.005488 

0.005*79 
S  005472 
:■  0  05  455 
5460 
0.005455 
0.005450 
0.005446 


31.6 

IS!. 

181. 

182.2 

182.5 

182. 

182.9 

183. 

183.; 

183. • 
183.6 


09 

01*5 
9925 
2054 
0460 
5151 
6809 
5962 
3020 
8309 
O90 


5.594 

181.93645 

5.5£~ 

181.50081 

5.531 

181.88961 

5.574 

132.21991 

5.568 

182.50597 

5.561 

182.75037 

5.555 

132.96746 

5.548 

183.15399 

5.542 

133.32957 

5.535 

183.48246 

5.529 

183.62027 

TWO-STATE   ENVIRONMENT: 


TEST  EXP.    DOWN 

PARAMETER  TIME/YEAR 

mu(ss)  (   HOURS   ) 


FAILURE  MEAN  MEAN  MEAN 

FREOUENCY       CYCLE   TIME  DOWN  TIME  UP  TIME 

(1/YEAR)         (   YEARS   )  (   HOURS   )         (  YEARS   ) 


18.338910 

0.039901 

0.005524 

181.02392 

26.172801 

0.030786 

0.00551O 

181.48788 

22.006692 

0.03O684 

0.005498 

181.87634 

23.840583 

0.030593 

0.005488 

182.20635 

25.574474 

0.030510 

0.005480 

182.49017 

27.508365 

0.030433 

0.005472 

182.73686 

29.342256 

0.035362 

0.005466 

182.95326 

31.176147 

0.005460 

183.14462 

33.910033 

0.93O232 

0.035455 

183.31  506 

34.843929 

0.030172 

0.005451 

183.46782 

36.57"?S29 

a  <f\a-\  i  ± 

9  ?65"6 

in;  ep»^«;-7 

WO-STATE  E 

NVIRONMENT: 

x  =  0.80 

5 

.594 

131 

02328 

5 

587 

131 

48724 

5 

.581 

181 

87571 

5 

.574 

182 

20571 

5 

.568 

1S2 

48953 

5 

.561 

182 

73623 

5 

.555 

182 

.95263 

5 

.548 

183 

14399 

5 

.542 

153 

.31443 

5 

.536 

133 

.46719 

*; 

<^*?o 

133 

c.o>  400 

TEST  EXP.    DOWN 

PARAMETER  TIME/YEAR 

mu(ss)  (   HOURS   ) 


FAILURE  MEAN  MEAN  MEAN 

FREOUENCY       CYCLE   TIME  DOWN  TIME  UP  TIME 

(1/YEAR)  (   YEARS   )  (   HOURS   )         (  YEARS  ) 


18.338910 

0.030939 

0.005531 

180.80253 

5.594 

180.80190 

28.172881 

0.030825 

0.005517 

181.26017 

5.587 

181.25953 

22.006592 

0.030"* 

0.005505 

181.64334 

5.581 

181.64270 

23  5*0583 

0.638633 

0.005495 

181.96835 

5.574 

181.96821 

25.574474 

0  050551 

0  005487 

182.24881 

5  568 

132.24818 

27.508355 

0  o;o*_5 

0.005480 

182.49216 

5.561 

182.49153 

29.342256 

0.030404 

0.005*73 

182.70564 

5.555 

182.70501 

31.  1  76147 

0.030338 

0.OO5468 

182.89443 

5.5*9 

132.89389 

33.010038 

0  0302"5 

0  005463 

183  06259 

5.542 

183.06196 

34  3*3929 

0 0  30215 

0  005458 

183.21332 

5.536 

183.21269 

36.677820 

0  0  30  153 

0.005454 

183.34920 

5.529 

183.34857 

TABLE   4. 


SOLUTION  OF   RELIABILITY   INDICES    FOR   DUAL  GENERATOR 
alpha  -   0.01  mulss)  VARIED 


ONE   STATE   ENVIRONMENT: 

TEST  EXP.  DOWN 

=  ^iM"Er         TIME      EAR 
-     JRS    ) 


FAILURE  MEAN  MEAN  MEAN 

FREQUENCY       CYCLE  TIME         DOWN  TIME         UD  "'ME 
(1/YEAR)         :    YEARS    .         I    HOURS   )         (  YEARS 


18.53891  8 

8.006318 

0.000058 

• 7364  152' 3 

5.529 

17364.45150 

23. '72861 

8.060516 

0.008057 

'74-3.78461 

5.516 

1  -4-;.7£2?5 

22  386692 

0. 000313 

0  :•  b  e  0  5  ~ 

•"555  "8154 

5.504 

1-565.78091 

23.849583 

0.90031  1 

0.000657 

'"644  25247 

5.491 

17644.26185 

25.6-4474 

0.000303 

0.000055 

'  "12.00162 

5.478 

17712.00099 

27.503355 

0.000308 

0.000056 

17771.06259 

5.466 

17771.06196 

29.342255 

0.000305 

0.06O056 

17823.01243 

5.453 

17823.01181 

31.176147 

0.000304 

0.000056 

17869.06217 

5.441 

17869.06154 

33.810038 

0.000303 

0.000056 

17910.16265 

5.429 

17910.16203 

34.843929 

0.000302 

0.000O56 

17947  07087 

5.416 

17947.07025 

36.577820 

0.000301 

0.000056 

17980.39697 

5.404 

17986.39635 

TWO-STATE   ENVIRONMENT: 


6.20 


TEST  EXP.   DOWN 

PARAMETER  TIME/YEAR 

mu^ss)  (   HOURS   ) 


FAILURE  MEAN  MEAN  MEAN 

FREOUENCY       CYCLE   TIME  DOWN  TIME  UP   TIME 

(1/rEAR)         (   rEARS  )         (  HOURS  )         (  YEARS  ) 


18.333910 

0.0OO319 

0.000058 

17354.35471 

5.529 

17354. 354yb 

20.172801 

0.00O316 

0.000057 

17463.53553 

5.516 

17463.53490 

22.006692 

0.006313 

0.686057 

17555.48565 

5.504 

17555.40502 

23.846583 

0.00031  1 

0.000057 

17633.77980 

5.491 

17633.7791  7 

25.674474 

0.000309 

0.000656 

17701.42755 

5.478 

17701.42692 

27.508365 

0.000308 

0.660856 

17760.40972 

5.466 

17760.40910 

29.342256 

0.660306 

6.666856 

17812.29138 

5.453 

17812.29075 

31.176147 

0.080305 

3.008056 

17858.231  12 

5.441 

17858.28650 

33  016038 

0 000303 

a  n«fifl56 

1  7fmO  3?Q»7 

5  49Q 

1  7ROC1  J?SA5 

34.643323 

0.  0OO  jO  L 

O.OtfOODO 

1 7936. 1  966  i 

5.416 

l/3JO.(-t?!3 

36.577820 

0.000301 

6.888856 

17969.47599 

5.404 

17969.47537 

TWO-STATE   ENVIRONMENT:  x   =   0  80 

TEST  EXP.   DOWN  FAILURE  MEAN  MEAN  MEAN 

PARAMETER         TIME/YEAR  FREOUENCY       CYCLE  TIME  DOWN  TIME  UP   TIME 

mu(ss)  (   HOURS   )  (1/YEAR)  (   YEARS   )  (  HOURS   )  (   YEARS   ) 


18  333910 
20.172801 
22.006692 
23.348583 
25.674474 
27.508365 
29.342256 
31.175147 
33.010638 
34.843929 
36.677820 


0.000322 
0.006319 
0.00O316 
0.000314 
0 060312 
0.00031  1 

0.0003O9 

0.0003O8 
0.000306 
0.000305 
0.800304 


0  000058 

0  600058 
0.000058 
0.000057 
0  OO0O57 
0.000057 
0.00O057 
0.660057 
0  080056 
0  00OO56 
0  000056 


17193 
1730O. 
17390. 
17467. 
17533. 
17591. 
1  7641 
17686 
17727. 
17763 
17795 


77563 
57904 
46383 
15900 
37266 
1  1878 
92696 
97983 
20288 
33592 
97393 


5.529 

17193.77500 

5.516 

173O0.57841 

5.504 

17390.46320 

5.491 

17467.15338 

5.479 

17533.37284 

5.466 

17  591.1  1316 

5.454 

17641.92634 

5.442 

17686.97921 

5.429 

17727.20226 

5.417 

17763.33530 

5.405 

17795.97331 

TABLE  5. 


SOLUTION   OF  RELIABILITY   INDICES   FOR   SINGLE   GENERATOR 
alpha  =  9.95  mul'ss'l  VARIED 


ONE  STATE   ENVIRONMENT: 

~E3~               EXP     DOWN  FAILURE  MEAN  MEAN  MEAN 

SARAME"^          TIME.  x  EAR  FREOUENCI  CYCLE  TIME  DOWN  TIME          UP   T'ME 

mu(ss                    HOURS  rEAR)  (    YEARS    i         (  HOURS   )          (  YEARS    ) 


20.772801 
22.06559: 
23.840583 

:e  S7+474 

2"  ;:-2355 
29.342256 

31.1761  47 
33.818638 
34.843929 
36.677820 


0.1 

0.1 

0 

0 

8.1 

0.1 
0.1 
0.1 

0.1 


53713 
53219 

s::5E 

51994 
51545 
51317 
51006 
50711 
56428 


S.82 

B  •:■: 

0.02 
0.02 

0.82: 

0.02: 

0.02: 

0.02: 

0.02 
0.02 
0.02 


7511 


7407 
7356 


7248 
7226 
7207 


36.25934 
36.34851 
36.42317 
36.48659 
36.541  1  4 
36.58355 
35.53014 
36.66692 
36.69968 
36.72904 
36.75550 


5.594 

36.255-C 

5.587 

36.34787 

5.581 

36.42253 

5.574 

36.48596 

5.568 

36.54651 

5.551 

35.55"?: 

c  ccc 

36.52951 

5.548 

36.56629 

5.542 

36.69904 

5.535 

36.72840 

5.529 

36.75487 

TWO-STATE  ENVIRONMENT:      x  =  0.20 

TEST      EXP.  DOWN  FAILURE      MEAN  MEAN  MEAN 

PARAMETER    TIME/YEAR  FREOUENCY   CYCLE  TIME    DOWN  TIME    UP  TIME 

mu(ss)      (  HOURS  )  (1/YEAR)    (  YEARS  )    (  POURS  )    (  YEARS  ) 


18.338910 
20.172801 
22.005692 
25.849583 
25.67^-74 
27.505365 
29.342256 
31.175147 
33.818038 
34  8*3929 
36  5~7320 


.154237 
.153679 
.153185 
.152741 
152335 
.151961 
.151612 
.151284 
.150974 
.  15*679 
.150395 


9.92 
8.02 
9.02 
0.02 
0.02 
8.02 
0.92 
0.82 

0.02 

0.02 


7573 
7585 
7449 
7481 

7261 
7225 
7294 
7267 
7242 
7221 
7281 


36.26730 
36.35644 
36.43108 
36.49449 
36.54902 
36.59642 
36.63799 
36.67476 
36.70751 
36  73686 
36.76332 


5.594 

36.266-66 

5.587 

36.35589 

5  581 

36.43944 

5.574 

36.49385 

5.568 

36.54e2S 

5.561 

C  =,=<£. 

36.53736 

5.548 

36.67412 

5.542 

36.79688 

5.535 

36.73623 

5.529 

36.76269 

TWO-STATE  ENVIRONMENT:     x  =  9.80 

TEST      EXP.  DOWN  FAILURE      MEAN        MEAN        MEAN 

PARAMETER    TIME/YEAR  FREOUENCY   CYCLE  TIME    DOWN  TIME    UP  TIME 

nu(ss)      (  HOURS  )  (1/YEAR)    (  YEARS  )    (  HOURS  )    (  YEARS  ) 


18.333919 
20.172801 
22.006692 
23.349583 
25  674474 
27.508265 
29.342256 
31.175147 
33.018038 
34.843929 
36.577820 


.153 
.153 
.152 
.152 
.151 
.151 
.151 
.158 
.150 
.150 


71  1 
158 
669 
228 
825 
454 
188 
783 
475 
182 
901 


0.027*79 
e. 827412 
0.027356 
8.827309 


0.0 

0.027 


:24 
?03 


0.827176 
8.027152 
0.027130 
0.8271  1  1 


36.39142 
36.48816 
36.55447 
36.61759 
36  67183 
36  71907 
36  76047 
36.79708 
36.82969 
36.85893 
35  S8528 


5.594 

36.39678 

5.587 

36.47953 

5.531 

36.55383 

5.574 

36.61695 

5.563 

36.57124 

5.561 

36.71843 

5.555 

36.75984 

5.548 

36.79645 

5.542 

36.82996 

5.536 

36  85829 

5.529 

36  88465 

TABLE   6. 


SOLUTION   OF   RELIABILITY   INDICES   FOR   DUAL  GENERATOR 
aloho  =  0.05  mu(ss)  VARIED 


TEST      I 

<P  2CWN 

FAILURE 

MEAN 

MEAN 

MEAN 

paramete= 

TIME/YEAR 

FREOUENCY 

CYCLE  TIME 

DOM 

TIME 

UP  TIME 

muisst 

HOURS  ) 

EAR)    ( 

YEARS  )    ( 

HOURS 

E-FS 

1  E. 332910 

o .  ■:•  0792  9 

0  1  434 

697.291  71 

s  c  "0 

657.29108 

20.1  72801 

0.067864 

-.' .  0  0 1 4-2  6 

701.46885 

5.5'.6 

761.46822 

22.006632 

0.007807 

9  001418 

704.98278 

5.504 

"04.98216 

23.  349583 

0.007756 

0.001  412 

707.97980 

5.491 

"07.97918 

25.674474 

9.007710 

e. 001407 

716.56612 

5.478 

716.56549 

2". 508365 

0.007668 

0.001403 

712.820"  " 

5.466 

"•;.  32009 

23. 3*2255 

9.00"€23 

0.001339 

714.80354 

5.453 

714.89291 

31.176147 

0.007S93 

0.00  1  396 

716.56092 

5.441 

"ic.56030 

33.019e5S 

0.06755S 

0.001393 

718.12922 

5.429 

718.12860 

34.843929 

0.807527 

0.001390 

719.53738 

5.416 

719.53676 

36.677820 

0.007497 

0.001387 

720.80873 

5.404 

720.8081  1 

TWO-STATE    ENVIRONMENT: 


0.20 


TEST  EXP.   DOWN  FAILURE  MEAN  MEAN  MEAN 

PARAMETER  TIME/YEAR  FREOUENCY       CYCLE   TIME  DOWN  TIME         UP   TIME 

mu(ss)  (   HOURS   )  (1/YEAR)  (   YEARS   )         (   HOURS   )  (   YEARS   ) 


18.338910 
20.172881 
22.006692 
23.840583 
25.674474 
27.508365 
29.342256 
31.176147 
33.010638 
3  4  2  433^9 
36.677820 


0.907929 
0.O07864 
0.007807 
6.867756 
6.867710 
6.607668 
0.667629 
8.667593 
8.007559 

0.007497 


001 434 
001426 
001 418 
001412 
081407 
901 403 
001399 
001395 
991393 
RBI  Mfl 


0.001387 


697.29223 

5.529 

697.29160 

701.46798 

5.516 

761.46645 

704.97989 

5.564 

704.97846 

787.97458 

5.491 

707.97387 

718.55942 

5.478 

716.55380 

712.81234 

5.466 

712.81221 

714.79463 

5.453 

714.79481 

716.551  1  1 

5.441 

716.550*3 

718.11862 

5.429 

718.1  1801 

7ig  526 1 1 

5  4  *  c 

728.79685 

5.404 

720.79623 

TWO-STATE   ENVIRONMENT: 


x   =   0.80 


TEST      EXP.  DOWN     FAILURE      MEAN        MEAN        MEAN 
PARAMETER    TIME/YEAR     FREOUENCY   CYCLE  TIME    DOWN  TIME    UP  TIME 
mu(ss)      (  HOURS  )      (1/YEAR)    (  YEARS  )    (  HOURS  )    (  YEARS  ) 


18.333910 
20.172881 
22.006692 
23.840583 
25.674*74 
27.508365 
29.342256 
31 . 1 76147 
33.016638 
34.843929 
36.677820 


0.067930 
0.007865 
0.007808 
0.807757 
6.007712 
8.007670 
0.067631 
0.007596 
0.007562 
0.867536 
0.067500 


0.601  434 
0.081426 
0.001419 
0.001413 
0.081403 
0.001403 
0.0013S9 
0.661396 
8.681393 
8.001390 
0.091388 


697.25629 
701.39499 
784.87674 
787.84560 
716.48969 
712.64429 
714  60986 
716.35217 
717.90728 
719.30385 
720.56503 


5  529 

697.25566 

5.516 

701.39436 

5.564 

704.87612 

5.491 

707.84597 

5.478 

716.40906 

5.466 

712.64367 

5.453 

714.60924 

5-44  1 

716.35155 

5.429 

717.90666 

5.416 

719.30323 

5.404 

720.56441 

TABLE   7. 


SOLUTION   OF   RELIABILITY   INDICES   FOR   SINGLE 
alpha   -0.10  mu(ss:  VARIED 


GENERATOR 


ONE    STATE    ENVIRONMENT 


TEST      E 

<P  DOWN 

FAILURE 

MEAN 

MEAN 

MEAN 

PARAMETER 

TIME/YEAR 

rRE3LESC 

CYCLE  ~IME 

DOWN  TIME 

UF  TIME 

ITOJ  { ss  1      ( 

HOURS 

(1/YEAR 

kears 

HOI  RS 

18.338916 

0.308036 

0.055'?" 

:S.'564g 

5.594 

18.15585 

20  '  72801 

0.367012 

0.054545 

13.19873 

5.55" 

'2 

22.806692 

0.366059 

6.054842 

12.23409 

5.581 

23.840583 

0.305199 

0.054752 

"8.26414 

5.574 

18.26550 

25.674474 

0.304413 

0.054675 

1S.2S998 

5.558 

IS. 22934 

27  50836S 

0.303686 

0. -354608 

18.31244 

5.561 

18.31180 

29  342256 

0.3O3O07 

0.054549 

18.332--4 

5.55: 

18.33151 

31 . i  751 47 

0.3O2368 

0.054497 

13.34956 

5  548 

18.34893 

33.010038 

0.301762 

0.054451 

18.36508 

5.542 

18.56445 

34.843929 

6.301  184 

0.054410 

18.37899 

5.535 

18.37836 

36.677820 

0.300631 

0.054373 

18.39152 

5.529 

18.39089 

TWO-STATE    ENVIRONMENT:  x    =    0.20 

TEST                EXP     DOWN  FAILURE                MEAN  MEAN  MEAN 

PARAMETER          TIME/YEAR  FREOUENCY       CYCLE   TIME  DOWN   TIME  UP   TIME 

MU(SS)              (   HOURS   )  (1/YEAR)         (   YEARS   )  (   HOURS   )          (  YEARS   ) 


18.338910 
20.172801 
22.006692 
23.840583 
25.674474 
27.508365 
29.342256 
31.176147 
35.010058 
i  *    a  4  loin 

36.677820 


308006 
306935 
305980 
505121 
504355 
505608 
502930 
302291 
301685 


i_\5l'05*5 


0.055O62 
0.054935 
0.054828 
0.054738 
0.054661 
0.054594 
0.054535 
0.0544E3 
0.054437 
a   £k  4305 

0.054359 


18.16120 
18.20343 
18.23879 
18.26883 
18.29467 
18.31712 
18.33682 
18.35424 
18.36975 
i  o  2325£ 

18.39619 


5.594 

8. 16056 

5.587 

8.20280 

5.581 

8.23816 

5.574 

8.26820 

5.568 

8.29403 

5.551 

8.31649 

5.555 

8.33619 

5.548 

3.35361 

5.542 

8  3£5"2 

ft  cie 

q  TB^cn 

TWO-STATE  ENVIRONMENT:  »  =  0.80 

TEST  EXP.  DOWN  FAILURE      MEAN        MEAN  MEAN 

PARAMETER    TIME/YEAR  FREOUENCY   CYCLE  TIME    DOWN  TIME    UP  TIME 

mu(ss)  (  HOURS  )  (1/YEAR)    (  YEARS  )    (  HOURS  )    (  YEARS  ) 


18.338910 
20.172801 

22.005692 
23.840585 
25.674^-1 
27.508565 
29.542256 
51.176147 
55.010058 
54.843929 
36.677820 


0.306759 
0.305695 
0.304751 
0.303899 
0.3031  19 
0.302398 
0.301725 
0.301091 
0  300490 
0.29991  7 
0.299367 


0.054839 
0.054713 
0.054608 
0.054519 
0.054442 
0.054376 
0.054518 
0.054267 
0.054221 
0.054181 
0  0541 44 


8.25505 
8.27714 
8.31238 
8.34252 
8.36807 
3.39045 
8.41009 
8.42746 
8.44292 
8.45679 
3  46929 


5.594 

8.23441 

5.587 

8.27650 

5.581 

8.31  175 

5.574 

8.34168 

5  558 

8.36744 

5.561 

8.38932 

5.555 

8.40946 

5.548 

8.42682 

5.542 

8.44229 

5.535 

8.45616 

5.529 

8.46866 

TABLE  8 


SOLUTION   OF  RELIABILITY   INDICES   "OS   DUAL   GENERATOR 
clpho  =   0.16  mu(ss)  VARIED 


ONE  STATE  ENVIRONMENT 

TEST      EXP.  DOWN 

rAILLRE 

MEAN 

MEAN 

MEAN 

=  ARAWE"  = 

TIME/YEAR 

FREOUENCY 

CrCLE  TIME 

row. 

TIME 

UP  TIME 

inu(ss)     ( 

HOURS  ' 

(1/YEAR)    ( 

YEARS  ) 

i 

rEARS 

18.338S  '  '• 

•3.03:534 

0.005712 

175.  9567 

5  =  ~5 

■~r      J5613 

28.1  "  2881 

0.031336 

9.065551 

176.63945 

5.516 

176.63832 

22.066892 

0.031117 

0.005554 

176.86587 

5.504 

176.86524 

23.840583 

?.e5??;3 

0.005532 

177.57855 

5.491 

1"". 56993 

25.674474 

0.033747 

0.805512 

178.17854 

5.478 

"8.17-52 

27.506555 

8.030585 

0. 005596 

178.76846 

5  466 

178.70783 

29.342256 

0.030A36 

0.005581 

179.174*2 

5.453 

179.17380 

31.176147 

6.030297 

0.805568 

".  79.587J4 

5.-14  1 

179.58672 

33.016938 

0.030166 

0.005557 

179.95580 

5.429 

175.95518 

34.843929 

0.030043 

0.005547 

180.28659 

5.416 

180.28597 

36.677820 

0.029925 

0.005538 

180.58521 

5.404 

180.58450 

TWO-STATE    ENVIRONMENT: 


0.20 


TEST  EXP.   DOWN 

PARAMETER         TIME/YEAR 
mu(ss)  (   HOURS   ) 


FAILURE  MEAN  MEAN  MEAN 

FREOUENCY        CYCLE  TIME  DOWN   TIME  UP  TIME 

(1/YEAR)  (  YEARS   )  (   HOUPS   )  (   YEARS   ) 


18.338910 

0.031582 

0.005712 

175.07216 

20.172801 

0.031333 

0.005680 

176.05443 

22.006692 

0.031  1  15 

0.005654 

176.33050 

23.840583 

0.030920 

0.805531 

177.58488 

25.574.17  a 

0.030744 

0.005612 

178.19261 

27.50e365 

0.036533 

0.005595 

178.72:31 

25.:-*22c6 

0.030434 

0.005531 

179.18808 

31.176147 

0. 030295 

0.005553 

179.60634 

-i  J'Pflio 

e.230'54 
0.030041 

0.005557 
8.805545 

■\  70  OCQ  T  5 

34.845929 

180.29581 

35.677820 

0  (IMO^l 

0.005537 

180.59335 

5.529 

175.07153 

5.516 

176.05380 

5.504 

176.87587 

5.491 

177.58425 

5.478 

178.19199 

5.466 

178.72168 

5.455 

175.18746 

5.44  1 

175.60022 

5  429 

1  "7Q  O  C  OC  ~t 

5  416 

180. 29920 

5.404 

ie0.5ST71 

TWO-STATE   ENVIRONMENT: 


0.80 


TEST      EXP.  DOWN      FAILURE      MEAN        MEAN        MEAN 
PARAMETER    TIME/YEAR     FREOUENCY   CYCLE  TIME    DOWN  TIME    UP  TIME 
mu(ss)      (  HOURS  )      (1/YEAR)    (  YEARS  )    (  HOURS  )    (  YEARS  ) 


'3.558910 
20.172301 
22.086692 
25.840583 
25.674174 
27.508365 
29.342256 
31.176147 
35.010658 
54.843929 
35  677820 


031 
031 
031 
6  3C- 
050 

>:■  5  ■:■ 

03O 
05 
838 
030 

629 


0.005704 
O.O05673 
0.005646 
0.005524 
0.005605 
0.005589 
0.005574 
0.905562 
0.005550 
0.605540 
0.005531 


75.3 
76.2 
77.1 
7"  8 
78.4 
78.9 
79.3 
79.8 
80.1 
80  4 
80.7 


0623 
8179 
0226 
9191 
0560 
3182 
9460 
0474 
7078 
9944 
9621 


c  5*70 

175.30568 

5.516 

176.281  16 

5.584 

177.10165 

5.491 

177.83129 

5.478 

178.40497 

5.466 

178.93120 

5.453 

179.59558 

5.44  1 

179.80412 

5.429 

180.17016 

5.416 

180.49835 

5.404 

180.79559 
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APPENDIX 

A  fortran  coded  program  for  the  solution  of  five  key  reliability 
measures  follows.  These  measures  can  be  solved  on  a  repeated  basis 
while  varying  a  selected  input  parameter.  Program  comments  provide 
information  regarding  data  input.  A  sample  data  input  file  for  the 
single  engine,  two-state  environment  (;■<  =  0.20}  is  provided  at  the  end 
of  the  program.  Equations  (1)  -  (7)  were  implemented  within  the  prograrr 


c 

C  •   •  CALCULATION  OF   RELIABILITY   INDICES   FOR         •   • 

C  •   «  EMERGENCY  AND   STAND-BY  POWER   SYSTEMS  •   • 

C 

C 

C  FUNCTION:   PROGRAM   WILL  PERFORM.  AS   DESIRED.   EITHER   A  ONE   TIME   CAL- 

C  CULATION   OF  RELIABILITY   INDICES   FOR   A   SYSTEM   BASED  ON   DATA 

C  AVAILABLE   FROM   THE   SYSTEM'S   STATE   TRANSITION   DIAGRAM  OR   WILL 

C  REPETITIVELY   COMPUTE  THESE   INDICES   WHILE  VARYING  A  SELECTED 

C  TRANSITION   RATE   OR   COEFFICIENT  OVER  A   SPECIFIED   RANGE. 

C 

C  DATA   INPUT  REQUIREMENTS: 

C 

C  A)   FIRST  DATA  ENTRY   IS   :  (INTEGER) 

C  B(ZERO)   =>   IF   ONE  TIME   COMPUTATION   DESIRED 

C  1  =>   IF  REPEAT  COMPUTATIONS   DESIRED 

C  NOTE(l):    IF   ZERO(0)   IS   ENTERED  THEN   DATA  INPUT   STEP(F)   MUST  BE 

C  SKIPPED. 

C  B)   NUMBER  OF   SYSTEM   STATES   (INTEGER) 

C)   NUMBER  OF  DISTINCT  TRANSITION   RATES  AND  TRANSITION   RATE 
C  CCE-r.C;ENTS    .INTEGER). 

NOTE.':):     «-   <<   MUST  COUNT    1.0   AS   DISTINCT   COEFFICIENT   >>   •-- 
NCTE.:      >~  THE  PROBABILISTIC   COMPLEMENT  OF  -   COEFFICIENT  IS  TO 
C  BE   USED.  "rHE   COMPLEMENT   NEED   NOT   BE   COUNTED   (IT  WILL  5E   INDICATED 

C  USING  A   MINUS   SIGN   IN   INPUT   STE~>  G) 

C  D!    TOTAL  NoMSE1    OF  TRANSITIONS   BETWEEN   STATES    'INTEGER! 

E     DISTINCT  TRANSITION   RATES    AND   COEFFICIENTS-   (i).(ii) 
INPUT   ORDER   nomEE^   (INTEGER) 

TRANSITION   rA7E  IN  UNITS   Cr    "       EAR  OR  UNITLESS  COEFFICIENT 
C  REAL) 

C  NC~E    I)     THE   r=ST  ENTRK    MUST  5E  THE  COEFFICIENT    1.9. 

C  F)   REPETITION   CONTROL 

C  '   INPUT  STEP(A)   IS    ".  THEN   ENTER:   (i).(ii).(iii) 

C  i)  THE   NUMBER   Or   THE   PARAMETER   (TRANSITION   RATE   OR   COEFFICIENT) 

C  TO   BE   VARIED   (INTEGER) 

C  ii)   INCREMENT   BY  WHICH   IT  IS   TO   BE  VARIED   (REAL) 

C  iii)   MAXIMUM  VALUE   FOR   PARAMETER   (REAL) 

C  NOTE(1):   MINIMUM  VALUE   OF  THE   PARAMETER   IS   ENTERED   AS   THE   VALUE   IN 

C  STEP(E). 

C  NOTE(2):    SKIP  THIS   STEP   IF   STEP(A)   INPUT   IS   ZERO(0). 

C  G)  TRANSITIONS:   (i).(ii).(iii).(iv).(v) 

C  i)   START   STATE   (INTEGER) 

C  ii)   END   STATE   (INTEGER) 

C  iii)   TRANSITION   RATE   NUMBER   (INTEGER)   -   CORRESPONDING  TO  NUMERICAL 

C  ORDER   RATES   WERE   ENTERED   EARLIER 

C  iv)  FIRST   COEFFICIENT   NUMBER   (INTEGER)   -   INDICATE   PROBABILISTIC 

C  COMPLEMENT  OF   COEFFICIENT  BY   PLACING  NEGATIVE   SIGN   IN   FRONT 

C  OF  THIS   NUMBER. 

C  V)   SECOND   COEFFICIENT   NUMBER   (INTEGER)   -    INDICATE   COMPLEMENT   IN 

C  SIMILAR   MANNER 

C  H)   NUMBER   OF   FAILURE   STATES   (INTEGER) 

C  0  LIST  OF  FAILURE  STATES  3r    STATE  NUMBER  (INTEGERS) 

C 

IMPLICIT  REAL«e(A-H.C— 2).  INTEGER(I-N) 

REAL»8  TRC(20).   B(A0.40),   P(40).   ZA(40,80). 
•TEMP(40),   BINV(40.40) 

INTFGEP  FS(*0),   FM(*0),   JSTPT(100).  JEND(100), 
>JTN(160),  vICN1(I66j.  JCN2(100) 

OPEN   (UNIT-1,   F!LE  =  'FELCUT.DAT',   STATUS-'NEW') 

CPEN   (UNIT-2.   FILE  =  'REL!N.DAT\   STATUS  =  'OLD') 

WRITE   (1.10) 
10   FORMAT   (//.1X,   'SOLUTION  OF  RELIABILITY  INDICES   FOR  './/) 

C 

C INITIALIZE  TRANSITION   RATE   MATRIX   (B-MATRIX) 

C 

DO   20   K-1.NS 
DO  23  L-1.NS 
B(K.L)   =   0  0 
20   CONTINUE 
C 

C   -   -    -      READ   IN   DATA 
C 

READ   (2.»)   IFLAG.NS.NTRC.NT 

READ  (2.»)  (M.TRC(I),I=1.NTRC) 

IF   (IFLAG.E0.1)  THEN 

READ   (2,«)   NVAR.H.PARMAX 

ENDIF 

READ  (2.«)  (JSTRT(I).JEND(I).JTN(I).JCN1(I).JCN2(I).I=1.NT) 

READ   (2.»)   NFAIL 

READ   (2,«)   (FS(I).  1-1. NFAIL) 
C 

C   -    -   -      PRINT  OUTPUT   HEADER 
C 

IF   (IFLAG.EO.0)  THEN 

WRITE   (1.23) 
23   FORMAT  (2X, 'EXPECTED   DOWN'. 5X. 'FAILURE'. 5X.3(5X, 'MEAN'. 6X). 
•/.2X/TIME   PER   YEAR'. 5X.'FREOUENCY'.6X. 'CYCLE  TIME'.6X. 

•  'DOWN  TIME'. 7X. 'UP   TIME',/.4X.'(   HOURS   )'.8X.'(  1/YEAR)',4X, 
•3X.'(   YEARS   )',3X.3X.'(   HOURS   )'.3X,3X.'(   YEARS   )'.3X. 

./.  80C  -■).//) 
ELSE 
WRITE   (1.27) 
27    FORMAT  (5X. 'TEST'. 7X, 'EXP.   DOWN'. 6X. 'FAILURE'. 3X. 

•  3(4X.'MEAN',5X)./.3X.'PARAMETER'.4X.'TIME/YEAR',5X, 

•  'FREQUENCY'. 3X. 'CYCLE  TIME'. 4X. 'DOWN  TIME'. 4X. 'UP  TIME'. 


•  /.16X.-(   HOURS   )'.6X. '(I/YEAR)-. 2X,2X.-(   YEARS   )\2X. 

•  2X.-(   HOURS   )-.2X.2X.-(  YEARS   )\2X./.80C-')./7) 
EN&lF 

C 

C FOPMULATION   OF  TRANSITION   RATE   MATRIX   (B-MATRIX) 

C 

30   DO   35   1=1. NT 

IF  (JCNI(l).LT.e)  THEN 
JCMPL   -    -JCN1(I) 
COEFF1    =    1.6   -   TRC(JCMPL) 
ELSE 

COEFF1    =   TRC(JCN1(I)) 
ENDIF 

IF   (JCN2(i).LT.0)  THEN 
JCMPL   =    -JCN2(I) 
COEFF2   -   '1.0  -   TRC(JCMPL) 
ELSE 
COEFF2   -  TRC(JCN2(I)) 
ENDIF 

B(JEND(I   ,JSTRT(I))   -  TRC(JTN(I))   .   COE~  1    -   CSl-'Z 
35   CONTINUE 
DC    40    I--1.NS 
8(1,.)  =  0.0 
DO  40  J»  i  .NS 

ecu    joto  *0 

B(l.      — B  I.I)  -    : 

4-0   CONTINUE 
C 

C  SAVE    ^CW    1    IN  TEMPCRART  ARRAY   ) 

C 

DO   50   1=1. NS 
TEMP(I)   =   B(1,l) 
B(l.l)  =    1.0 
50   CONTINUE 
C 

C   -    -   -      INVERT   8-MATRIX 
C 

NP    =    NS    +    NS 

CALL   MATINV(NS.NP.B.BINV.ZA) 
C 

C   -    -    -      COMPUTATION   OF   STEADY   STATE   PROBABILITY  VECTOR   (P-VECTOR) 
C 

DO   60   1=1. NS 
P(l)    =    SINV(I.I) 
60   CONTINUE 
C 

C   -    -   -     RESTORE  ORIGINAL  TRANSITION   RATE   MATRIX   (B-MATRIX) 
C 

DO  70   1=1  .NS 
B(1.l)   =   TEMP(I) 
70   CONTINUE 
C 

C 

C 

C       •••   COMPUTATION   OF   RELIABILITY  INDICES   «•• 

C 

C FORMULATION  OF   FAILURE   MODE  VECTOR   (FM-VECTOR) 

C 

YTH   =    8766.9 
DO   80   1=1. NS 
FM(I)   =   0 
80   CONTINUE 

DO   99   1=  1.NFAIL 
K    =    FS(I) 
FM(K)   =    1 
98   CONTINUE 
C 

C   -    -   -<<<   EXPECTED   DOWN   TIME   PER   YEAR   (HOURS)   >>> 
C 

PF   -   0.0 

DO    100   1=1. NS 

IF   (FM(I).EO.I)  THEN 

PF   -   PF   +   P(l) 
ENDIF 
100   CONTINUE 
C 

C <<<   FAILURE   FREQUENCY   OF   SYSTEM   (FF)   >>> 

C 

FF   =   O.0 

DO    120   1=1. NS 

IF   (FM(I).EO.O)   THEN 

SUM    =    0.0 

DO   1  10  J=  1  .NS 

IF   (FM(J).EO.I)   THEN 
SUM   =   SUM    +    B(J.I) 
ENDIF 
110       CONTINUE 

FF   =   FF  +   P(I)»SUM 
ENDIF 
120   CONTINUE 
C 

C   -    -    -<<<   MEAN   CYCLE  TIME   (TF)  /  MEAN   DOWN   TIME   (DF)   /  MEAN   UP   TIME   (DU)   >>> 
C 

TF   =    1.0/FF 
DF   =   PF/FF 


DU   -   TF   -   DF 
C 

C 

C 

C     «•   OUTPUT   .. 

C 

PF   =    PF.YTH 
CCC  TF   =   TF.YTH 

DF   =    DF.YTH 
CCC  DU   =   DU'YTH 

IF  (IFLAG.E0.1)   THEN 

WRITE   (1.139)  TRC(NVAR),PF.FF.TF.DF.DU 

139  FORMAT  (1X.3F13.6.F13.5.F13.3.F13.5) 
IF  (TRC(NVAR).LT.PARMAX)  THEN 

TRC!NVAR)   =   TRC(NVAR)   -i-   H 

GOTO   39 
ENOIF 
ELSE 

WRITE   (1.148)   PF.FF.TF.DF.DU 

140  FORMAT  (1X.2F16.6.F  15.5.F13.3.F13.5) 
ENDIF 

STOP 
END 


SUBROUTINE   MA7INV<N.NP.A.AI.~A) 
C 

C        NOTE  N   =   DIMENSiON  OF  MATRIX   TO  BE  INVERTED   (   LESS  THAN  OR 

C  ECUAL  T0   DIMEvJSiCN   OF   MATRIX   AS   SPECIFIED   IN  THE 

C  DIMENSION   STATEMENT  OF  THE   CALLING  PROGRAM   ) 

C 

C  NP   =   2   •   N 

C 

C  A   =    MATRIX   TO  BE   INVERTED 

C 

C  Al   =    INVERTED   MATRIX 

C 

C  ZA    =    DUMMY   ARRAY 

C 

REAL«8    A(40,40).ZA(40,80).AI(40,40) 

REAL<8  ALFA, BETA 

DO   26   1=1, N 

DO   24   J=1.N 

ZA(I.J)   =   A(I.J) 

NJ   =   N   +   J 
24  ZA(I.NJ)    =    0.0 

Nl    =    N    +    I 
26  ZA(I.NI)    =    1.0 

DO  111    IR  -   I.N 

IF  (DABS(ZA(IR.IR)).LT.1.E-20)THEN 

DO   IROW  =  1  ,N 

IF(DABS(ZA(IROW.IR)).GE.1.E-20)THEN 

DO   IC   =    1.NP 

ZA(lR.IC)   =   ZA(IR.IC)   +   ZA(IROW.IC) 

cur,    r\r\ 

GO  TO    111 

END   IF 

END   DO 

GO  TO  35 

END   IF 
1  1  1  CONTINUE 

C 

C     «...   SET   MAIN   DIAGONAL  ELEMENT  TO   UNITY  AFTER 
C     ••••   TESTING  FOR   ZERO  VALUE 
C 

DO    12   1=1.  N 

ALPHA   =    ZA(I.I) 
C  IF(DABS(ALPHA).LT.1.E-20)GO  TO   35 

DO   5  J=1.NP 
5  ZA(I.J)   =    ZA(I.J)/ALFHA 

C 

C     <•••   SET  THE   ELEMENTS   CF  THE   ITH   COLUMN  TO   ZERO 
C 

DO    11    K=  1  ,N 

IF   (K  — 1)8. 1  1.8 
8  BETA   =    ZA(K.I) 

DO    10  J  =  I.NP 
10  ZA(K.J)   =    ZA(K.J)   -   BETA  •    ZA(I.J) 

1  1  CONTINUE 

12  CONTINUE 

DO   33  J=1.N 

JN    =  J   +   N 

DO   33   1=1. N 
33  AI(I.J)   =    ZA(I.JN) 

RETURN 
C 

C     ....   PRINT   ERROR   MESSAGE   IF  METHOD   FAILS 
35  PRINT   .,-....>        MAIN   DIAGONAL  ELEMENT   IS   ZERO   ••••>■ 

RETURN 

END 
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